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ABSTRACT 



Tables of mean opacities are often used to compute the transfer of radiation 
in a variety of astrophysical simulatio ns from stellar e volut ion models to proto- 



planetary disks. Often tables, such as lFerguson et al.l (120051 ). are computed with 



a predetermined set of physical assumptions that may or may not be valid for 
a specific application. This paper explores the effects of several assumptions of 
grain physics on the Rosseland mean opacity in an oxygen rich environment. We 
find that changing the distribution of grain sizes, either the power-law exponent or 
the shape of the distribution, has a marginal effect on the total mean opacity. We 
also explore the difference in the mean opacity between solid homogenous grains 
and grains that are porous or conglomorations of several species. Changing the 
amount of grain opacity included in the mean by assuming a grain-to-gas ratio 
significantly affects the mean opacity, but in a predictable way. 

Subject headings: dust, extinction — equation of state — methods: numerical - 
astronomical data bases: miscellaneous 



1. Introduction 

The opacity of gaseous environments is a topic of interest in many areas of astrophysics 
including stellar atmospheres, circumstellar and protoplanetary disks, and the interstellar 
medium. Often, at least in the case of stellar evolution calculations, opacity tables are 
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needed at many temperature, density and chemical composition points. For high temperature 
opacities modellers often use either the OP AL opacities ( Iglesia s & Rogers (1991, 1993, 
1996), Rogers k, Ig lesias (1992a, 1992b) and iRogers et al.l (119961 )) or the Opacity Project 
( ISeaton et al.l 11994 OP hereafter) tables. We have produced low temperature opacities at 
Wichita State University ([Ferguson et al.ll2005l F05 hereafter) since 1975. These tables differ 
from their higher temperature counterparts by including the effects of molecules and dust in 
the equation of state (EOS hereafter) and the opacity. 

Typical opacity tables provided by our group include over 1600 temperature and density 
points with 155 chemical compositions needed for a complete set. Given the large number of 
input variables, atomic and molecular physics, grain parameters for grain physics, etc., it is 
not computationally possible to make a complete set of opacity tables for every astrophyisi- 
cally significant situation. Difficult decisions must be made as to how to choose important 
physical parameters. For example, in the opacity tables of F05 grains are assu med to have 
a size distribution like the one described in a classic model (IMathis et al.l 1 19771 . MRN here- 
after) based upon grain sizes in the interstellar medium (ISM); however this choice may not 
be appropriate for different astrophysical environments. This paper discusses some of the 
choices made in F05 and how those choices affect the total Rosseland mean opacity. 

The calculations discussed here are taken from the same code used to complete the 
opacity tables described in F05, which includes a summary of the equation of state and 
opacities. For consistency we begin with the total Rosseland mean opacity, which is defined 

as 



j_ __ Jo k x 9 dT ^ /,X 

Kr Jo -ar dX 

where k\ is the monochromatic opacity and 8B\/dT is the derivative of the Planck func- 
tion with respect to temperature. The monochromatic opacity includes the effects of all 
contributors: 



K A — f^cont + K-at + ^mol + ^gr (2) 

and includes continuous, atomic, molecular and grain sources of opacity. Note that the 
Rosseland mean is a harmonic mean, one in which individual contributors cannot be added 
after the mean is taken. This is particularly true of the molecular opacities which may have 
a large number of individual lines or bands. If it is desired that an individual source of 
opacity be modified, it is usually the case that the whole integral must be recomputed. A 
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detailed description of these sources is given in F05. We focus in this paper on the details 
of the grain contribution to the monochromatic opacity. 

The opacity due to a particular dust species can be computed from 

K gr p = n i( a ) Qcxt(a,i,X)a 2 da (3) 

i J a 

where rii(a) is the number of dust particles (cm~ 3 ) of species i of size a and Q e ^t(a, i, A) is the 
total extinction (absorption plus scattering) efficiency of the particle. The size distribution 
rii(a) depends upon both the number abundance of species i and on the size distribution of 
the dust particles. Ordinarily in our opacity tables the size distribution is taken from MRN, 
the classic distribution of dust sizes in the ISM given by 



rii(a) = ka q 



(4) 



where k is a normalization constant of the distribution and q is the exponent of the power 
law. In F05 the value of q is taken to be the classical MRN value of -3.5. 

Additional work by many other groups have resulted in other "favored " size distribu- 
tions, however none have fully superseded the work of MRN (jClayton et al.l (120031 )). Other 
distributions such as Kim et al. (1995, KHM hereafter) and Weingarnter & Draine (2001, 
WD hereafter) use power-laws as their basis, but include transitions for large g rains so that 
the di strib ution behave s smoo thly. For carbon grains WD, using the work of iDraine &: Li 
( 120011 ) and lLi fc Draind (120011 ) go a step further and include multiple distributions for differ- 
ent sizes of grains. While this is not an exhaustive list of grain size distributions used in the 
astrophysical community today, it is representative of the types of grain models available. 
Each of the distributions listed have been inserted into our code for comparison with the 
MRN distribution. 

In Eq. 3, the grain extinction efficiency is computed with Mie theory and assumes 
that grains are solid. This may not be the case in all, or any, astrophysical environments. 
Many researcher s argue that interstellar g r ains are likely t o be porous or " f luffy" in nature 
(see for example iMathis fc Whiffenl Jl989h . iMathisI Jl996h . Ichiar fc Tielensl (120061)1 Havin g 
grains whic h are p art v acuum can significan t ly affe ct their opacity (see IWolff et al.l (11994 ) , 
Wolff et al.l (119981 ) and IVoshchinnikov et al.l (120051 ) ) . However, in astrophysical grains the 
amount of porosity, generally d efined as the volume fraction of vacuum, is not well con- 
strained; IChiar fc Tielensl (120061 ) interpret observ ations to const rain the porosity of dust in 
the ISM to be in the 25% to 50% range whereas IMathisI (119961 ) indicates that dust mayby 
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up to 80% vacuum. 

Related to the porosity of grains is whether or not grains are composed of aggregates or 
are homogeneous. It is likely that astrophysical grains are not homogeneous as is assumed in 
the opacity tables of F05, but rather are an aggregate of two, three or mor e different species . 



The effect of heterogeneous grains on the mean opacity using the utilities of lOssenkopj (If 9911 ) 
will be discussed below. 

The opacity tables of F05 assume that gas and dust are in equilibrium, but this may 
not be the case in real environments due to dynamical effects such as winds or gravitational 
settling of grains. One simple way to mimic the effects of grains either settling or being blown 
out of an environment is by changing the grain-to-gas ratio by assuming certain percentages 
of grain material have been removed from the gas. 

This work is an exploration of the effects of grain physics on the total Rosseland mean 
opacity of a gas at coo l temperatures (be l ow ~ 3200 K). Opacities calculated here are for 



solar abundances, from iGrevesse fc Noelsl (119931 ) scaled to X=0.7 and Z=0.02 (where X is 



the mass fraction of hydrogen and Z is the mass fraction of metals), for 2.70 < log T < 3.30 
and constant logi? = —3.0 (see F05). The value of \ogR = —3.0 best models the atmospere 
of an AGB-type star and was arbitraily chosen here as a convenient way to demonstrate 
grain physics effects. All comparisons are done in an oxygen rich chemistry and include only 
silicate grains, future work will emphasize carbon grains. We will discuss five important 
aspects of grains physics in the following sections: (1) the size distribution exponent (see 
Eq. 4), (2) the size distribution itself, (3) grain porosity, (4) aggregate grains, (5) and the 
relative grain-to-gas ratio. 



2. Grain Physics 

In gas at high temperature, above ~ 5000 K, neutral and ionized atoms (including lines) 
are the dominant opacity sources. As the temperature decreases, molecular effects become 
significant and begin to dominate the opacity (see F05). Below ~1600 K at logi? = —3.0 
grains begin to appear in the equation of state and quickly become the strongest source of 
opacity. For purposes of illustration, a plot similar to one from F05 is given in Figure 1 
with the regions of dominant opacity marked. The features below logT = 3.2 are caused by 
various grain species appearing or disappearing from the equation of state. At logT = 3.2, 
solid AI2O3 appears in the EOS and causes a sudden rise in total opacity as the temperature 
decreases. Silicate grains begin to appear at log T ~ 3.08. The opacity has another bump 
at log T = 3.02 as solid Fe appears and becomes the dominant source of opacity. At 
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temperatures below log T = 3.0 iron and silicates are the predominate opacity sources 
(F05); no additional significant species appear or disappear and the mean opacity is fairly 
featureless as it gradually lessens as the temperature decreases. 

The opacity tables computed in F05 assumed spherical, solid, homogeneous grains with 
sizes that obey the MRN size distribution and with absorption and scattering cross-sections 
that can be simulated with Mie theory. We further assume in F05 that the grains form 
in equilibrium with and remain with the gas they form in. We will relax a few of these 
assumptions in the following sections to explore how the Rosseland mean opacity is affected. 

2.1. Size Distribution Exponent 

Equation 4 is a normalized power law distribution. As the value of the exponent, q, is 
varied - as the power law steepens or flattens - the size distribution will vary. As q increases 
the function will become steeper, that is, there are more small grains than large grains. 
The distribution is normalized for each grain species so that the volume of grain material is 
constant for the different values of q. The normalizations and summations are performed in 
terms of log(a), where a is the grain size. 

As an illustration we include a plot of the "guts" of Eq. 3 in Figure 2. The figure shows 
how Q(a)n(a)a 2 behaves for a single grain species (AI2O3) and a single wavelength (1 /im). 
As the exponent of the size distribution increases from -2.5 to -5.5, the size distribution, 
n(a), steepens, but Q(a)n(a)a 2 flattens out. This flattening is due to the combination of the 
a 2 factor and the extinction efficiency Q increasing with size (as the grains are smaller than 
the Rayleigh limit) for the species shown. 

What happens to the mean opacity if the value of the power law exponent, q, is varied 
from —2.5 to —5.5? Figure 3 shows that the mean opacity in the dust region (logT < 3.2) 
increases as the power law flattens out. Considering Fig. 2 this result is not surprising. As 
the grain size distribution, q, flattens out the extinction efficiency, which is a function of 
grain size but not distribution, rises for larger grains. Thus a flatter grain size distribution 
contains more large grains and those larger grains contribute more to the total opacity. 

The overall result for the mean opacity is that a flatter grain size distribution results 
in more opacity by as much as a factor of three when the power law exponent changes from 
—4.5 to —2.5. It is easy to see the difference changing the exponent, q, to the mean opacity 
with the lower panel of Fig. 3. In the lower panel the logrithm of the ratio of the mean 
opacity computed with a certain q value compared with the baseline value of —3.5 is shown. 
This type of plot easily allows for the change in the mean opacity to be seen. 
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Fig. 1. — Logarithm of the total Rosseland mean opacity as a function of temperature for 
solar abundances and logi? = —3.0. Regions where different components dominate the 
opacity are indicated. Low temperature features are discussed in the text. 
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Fig. 2. — Plot of Q(a)n(a)a 2 from Eq. 3 for Al 2 3 at a single wavelength (1 /im). The 
values shown for the power-law exponent are -2.5 (dotted), -3.5 (solid), -4.5 (dashed), and 
-5.5 (dash-dot). 
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Fig. 3. — Comparison of opacities calculated with different values of the MRN power law 
exponent, q, describing the power law size distribution of grains. In the upper panel values of 
q are, from top to bottom, —2.5 (dotted), —3.5 (solid), —4.5 (dashed), and —5.5 (dash-dot). 
In the lower panel is the differences in the logarithm of the mean opacities for the various test 
cases and the baseline (q = —3.5) are shown. The dotted line is the logarithm of the ratio 
of the q = —2.5 case and the baseline, the dashed line is q = —4.5 case, and the dash-dot 
line represents the q = —5.5 case. The thin solid line is included to help guide the eye and 
is the value of the logarithm of one. 
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There are several features seen in the lower panel of Fig. 3 that deserve discussion 
such as the stair step behavior of the comparison with the q = —5.5 computation with the 
baseline. We include in our computations several species of grains in our total mean opacity 
with their different opatical constants and extinction efficiencies. As different species appear 
or disappear from the EOS the mean opacity is affected in a variety of ways based upon the 
dominant species. More particularly, there is a large dip in the comparison with the q = —5.5 
mean opacity between 3.04 < logT < 3.08. This feature, and the much smaller one seen in 
the q = —4.5 computation, is due to the differing effects on changing the size distribution 
for different species. In the lower panel of Fig. 3, at logT = 3.08 the difference between 
the baseline and the q = —5.5 computation begins to move downwards as the Mg-silicates 
begin to dominant the opacity. The difference then moves upwards at logT = 3.04 as solid 
Fe becomes the dominant source of opacity (see F05). Changing the size distibution in the 
same way for every species in our EOS has a different effect on that species' contribution to 
the total mean opacity. 

Notice also that the comparisons with the q = —2.5 and q = —4.5 are not symmetrical 
about the reference line 0.0. This shows that changing the value of the exponent does not 
lead to a linear, or clearly defined change in the mean opacity This fact is not surprising by 
noting in Fig. 2 that the extinction efficiency, a function of grain size, rises with a, but not 
in a "linear" way. 

It would be convienient to have a scaling factor to convert the baseline opacities to a 
set of mean opacities with a different grain size distribution and this may be possible away 
from transistion zones. However, at those temperatures where the dominate grain species 
are changing a simple scaling factor is not possible. 



The size distributions discussed in Section 1 have been added to our opacity code and 
to be complete we include a short summary of the KMH and WD grain size distributions. 
The size distribution from KMH is given by 



where a& is the transition size for the exponential decay. The values chosen for our compu- 
tations are discussed below. 

A more sophisticated grain distribution, which includes a smoother transition from the 
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Grain Size Distributions 




(5) 
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power-law to the expontial decay is given by WD for silicate dust and has the form 



where k is a normalization constant, q is the power-law exponent (originally a in WD), a t 
and a c provide for the cutoff and shape of the distribution, and the power-law exponent 
represented by q. The curvature of the function (changed by modifying j3) is given by the 
term 



Similar distributions are assumed by WD for carbon or graphitic grains which include a 
log-normal distribution for the grain sizes to account for the very small grains thought emit 
the "Unidentified Infrared" (UIR) bands seen in the diffuse ISM. The comparisons in the 
present paper focus on oxygen-rich mixtures where carbon grains do not exist under these 
conditions. 

For our comparisons we use the values given by WD for silicate dust in the ISM, that 
is at =0.17 /im and a c =0.1 \xm with a value of /?=0.3 and a power-law exponent of -2.1 for 
all grain types. The computations using the KMH distribution use their stated parameters 
as well. That is, a& is taken to be 0.14 /im and a value of -3.06 for q is used for silicate dust. 

Figure 4 shows the comparsion between Rosseland mean opacities computed with grain 
size distributions from MRN, WD and KMH. In the figure the mean opacites are shown in 
the upper panel and the differences in the logarithm of the opacities is shown in the lower 
panel. The results show that the mean opacities are somewhat similar, with MRN having 
the lowest values, WD having higher values, and KMH having the highest values of mean 
opacity. It is interesting to note that KMH, with a very simple exponential transition for 
large grains, tend to tail upwards at the lowest temperatures. At logT = 2.7 (500 K), silicate 
grains are the dominant source of grain opacity, with solid iron grains falling in abundance 
in favor of FeS (see Fig. 1 of F05). 

There is also a strange series of bumps in the lower panel of Fig. 4 at intermediate 
temperatures, 3.00 < logT < 3.15. As discussed in the previous section, these temperatures 
are a transition region where the strongest opacity source is passing from minor grain species 
to the silicates to solid Fe. Different grain species have different contributions to the opacity 
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based upon which form of the size distribution is used due to different optical constants and 
extinction efficiencies. The feature seen in the lower panel of Fig. 4, centered at logT = 3.08, 
in the computations with KMH compared with the baseline is much more pronounced than 
the computations with WD. The difference between the KMH and WD grain distributions 
is how the size transition to large grains is treated. The minor grain species are much more 
sensitve to the KHM distribution than the solid Fe grains are, since when solid Fe dominates 
the opacity both the KMH and WD computations appear to be a simple scaling factor 
greater in the mean opacity. 

As discussed in the previous section, there does appear to be a simple scaling factor 
that might apply to the grain distibutions, but only for a limited range in temperature, 
2.80 < logT < 3.02. At other temperatures, where the strongest opacity source is changing 
from one grain species to another, no such simple scaling factor is evident. 

Overall the effect on the mean opacity is fairly significant with the KMH calculation 
differing from the baseline (MRN) by 0.27 dex at logT = 2.9 (800 K), or nearly a factor 
of 3. As Fig. 4 shows the difference from the baseline is not constant with temperature, 
especially in the regions where the significant contributors to the opacity are from different 
grain species (with different extinction efficiences) . For example around logT ~ 3.15 the 
differences between the computations are smaller than at cooler temperatures. This implies 
that arbitrarily adding a factor of 3 to the mean opacity tables of F05 to mimic the effects of 
different size distributions is not a good approximation. This is particularly true when com- 
bined with the mean opacity changes with differing size distribution exponents as described 
above. 



2.3. Grain Porosity 



Calculations of grain efficiencies due to grain porosity are taken from lBohren fe Huffman 

( 119831 ). To initially probe the effect of porosity on the mean opacity we assume that the 



vacuum inclusions are spherical and we use either the Maxwell-Garnett or the Bruggeman 
average dielectric constant for the grain/vacuum composite. Maxwell-Garnett is more ap- 
propriate for a mixture where the matrix and inclusions are well defined. The Bruggeman 
approximation is more symme tric and valid for a t wo- co mponent mixture with matrix and 
inclusions indistingishable. As lBohren fe Huffman! (119831 ) discuss both Maxwell-Garnett and 
Bruggeman are equally valid, but either one or the other will be more valid in a particular 
circumstance. 



In our computations we find that certain grain species misbehave when using the Brugge- 
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Fig. 4. — Comparison of the mean opacity due to differing grain size distrbution models. In 
the upper panel, the solid line is the mean opacity when using the MRN distribution, WD 
is represented by the dotted line, and the model of KMH is shown by the dashed line. In 
the lower panel the differences are shown as in Fig. 3, with the thin solid line at the 0.0 
axis shows no change to the opacity, the dotted line compares WD with the baseline, and the 
dashed line compares KMH with the baseline. See the text for detailed discussion. 
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man approximation. Most particularly troubling is solid Fe. At certain wavelengths the k 
portion of the optical constant will be unphysically negative, causing the mean opacity to 
increase discontinuously by large factors. To test our metho d we have utilized the publi- 
cally available average dielectric routines of lOssenkopil (Il99ll ) and found that program fails 
to compute porous Fe grains without either blowing up with a fatal error, or producing 
negative parts of the optical constants. We have decided to remove Fe from the porosity 
computation and let it remain solid and not porous. We assume that the other grain species 
are valid under the Bruggeman assumption, the EOS is mostly dominated by MgSi0 3 and 
Mg 2 Si0 4 at 1000 K and below (see F05, Fig. 1). 

Results for a range of grain porosities, denoted by the volume fraction /, are shown in 
Figure 5 with the mean opacities shown in the upper panel. The baseline value is /=0.0 
(solid grains) and the tested values include / = 0.1, 0.25, 0.5, and 0.8. Recent work by 
Chiar fe Tielend (120061 ) report porosities of grains in t he ISM as being between 25% and 
50% vacuum by volume, in line with previous studies by I Wolff et al.l (119931 ) . The lower panel 
shows the opacity differences compared with the baseline of / = 0. 

As the porosity of grains is increased, the the total mean opacity is mostly unaffected, 
except for the largest / value. This is somewhat intuitive in that the mean opacity is 
already dominated by small particles, making those small particles slightly larger, does not 
change the mean opacity greatly. At the largest value of porosity, / = 0.8 the mean opacity 
does change by as much as 50% at the lowest temperature shown. Note at "intermediate" 
temperatures, 3.05 < logT < 3.13, the lower panel of Fig. 5 clearly shows that the mean 
opacity decreases with increasing porosity. The materials responsible for most of the mean 
opacity at these temperatures are SiC>2, calcium-silicates, and MgAl 2 04 for the most part 
(see F05, Fig. 1). As these types of grains become larger and more porous that they also 
become more transparent, unlike the silicates at cooler temperatures, where the mean opacity 
slightly increases with porosity. 



2.4. Aggregrate Grains 



To explore the use of aggregate grains and their effect s on the Rosseland mean o pacity 
we compute the average dielectric function as outlined by iBohren fe Huffman! (119831 ) using 
effective medium theory (EMT). More specifically we assume Maxwell- Garnett theory for 
ag gregates which is valid as long as the inclusions obey the Rayleigh limit as illustrated 
by IWolff et al.l (119981 ) . We use this simple EMT approximation rather than a more formal 
solution such as a discrete-dipole approximation (DDA) due to limits on computational time. 
In addition, many of the features of a DDA are not needed for mean opacities computations, 
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Fig. 5. — Comparison of varying the porosity of grains. In the upper panel, the mean 
opacities are shown. The baseline calculation is the solid line with dust porosity, /=0.0, the 
dotted line is /=0.1, the dashed line is /=0.35, the dash-dot line is /=0.5, and the dash-dot- 
dot line is /=0.8. In the lower panel the differences are shown as in Fig. 3, with the porosity 
/=0.1 compared with the baseline as dotted, /=0.25 as dashed, /=0.5 as dash-dot, and /=0.8 
as the dash-dot-dot line. 
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such as the scattering and polarizaton phase functions. 

To compare with the baseline runs as outlined above we assume an aggregate grain 
that is similar in composition to the silicate species computed in our EOS at logT = 2.9 
(800 K) and with ratios similar to silicates found in t he inclusions of a typical meterorite 



such as Semarkona as discussed in lHewins et al.l (119961 ). While our EOS did not identically 



produce the Semarkona results, qualitatively it was fairly close; we compared our EOS with 
Semarkona and choose a median value. The ratios of enstatite, forsterite, fayalite, and 
silicon dioxide were assumed to be 42%, 36%, 12% , and 10% r espect ively. We computed the 



average dielectric functions using the routines of lOssenkopil ( 119911 ) with the stated ratios. 
Solid iron grains were kept separate from the aggregate as is found in meteoritic inclusions. 
The mean opacity was then computed at a single temperature and density (log R) point and 
compared with the baseline. Repeating this numerical experiment at lower temperatures did 
not dramatically change the results 

The comparisons show that aggregate grains have a lower mean opacity, on the order of 
0.05 dex (12%) below the baseline, which is a smaller change shown in Fig. 3 for changes in 
the grain size exponent q from -3.5 to -4.5. The amount of change in the mean opacity for 
aggregate grains from the baseline is very modest and near the estimated range of error in the 
grain thermodynamic and optical data of 0.02 dex (5%) based upon comparisons with other 
opacity databases (see F05). This result implies that including or not including aggregate 
grains in the computation of the mean opacity is not an important decision when computing 
mean opacity tables. 



2.5. Grain-to-Gas Ratio 

The opacities of F05 were not intended to be used in non-equilibrium conditions. How- 
ever, we can mimic such environments such as conditions of gravitational "rain-out" . If the 
gas/dust environment is under a gravitational influence with insignificant radiation pressure, 
then grains that should appear in the equaton of state will settle down into the gravitational 
potential and "disappear" from view. This will leave behind the apperance of a grain de- 
pleted, or a gas rich, environment. This is best simulated by letting the grains appear in the 
equation of state, but not counting their opacities. Figure 6 shows such simulations. The 
baseline amount from F05 includes 100% of the grain opacity for each species. The upper 
panel of Fig 6. shows the effects of including 50%, 10% and 1% of the total grain opacity. 

The effect on the mean opacity is roughly equivalent to the amount of grains removed 
from the grain opacity term. That is, if the grain opacity is cut in half, then the total mean 
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Fig. 6. — Comparisons of total Rosseland mean opacity with varying amounts of grain 
opacity included. In the upper panel, the solid line is the baseline calcuation (100% of the 
grain opacity included), the dotted line includes only 50%, the dashed line 10% and the 
dash-dot line only includes 1% of the grain portion of the opacity. In the lower panel the 
differences with the baseline are shown as in Fig. 3 with the dotted line comparing the 50% 
ratio to the baseline, the dashed line comparing the 10% ratio, and the dash-dot line the 1% 
ratio. 
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opacity is also roughly half of the baseline value because the grain opacity usually completely 
dominates the total opacity whenever grains are present. This is easy to see in the lower 
panel of Fig. 6. Note that this effect is not true for all temperatures however. At higher 
temperatures, the differences in the opacites (between the baseline and the test cases) are not 
smooth. Just below the grain appearance temperature at logT ~ 3.2 there is still a modest 
molecular opacity present. If a small amount of grain opacity, say 1%, is added to the (now) 
larger molecular opacity, the roughly linear scaling does not apply. As the molecular opacity 
drops at lower temperature, there is a roughly linear relationship between the amount of 
dust included in the monochromatic opacity and the total mean opacity. This result means 
that to a reasonable approximation the opacity tables of F05 could, in principle, be modifed 
to account for the effect of a grainless gas by multiplying by the desired factor, but only at 
the coolest temperatures. It would not be advisable near the transitions between molecules 
and grains dominating the opacity. 



3. Discussion 

This work implies that caution should be taken when requesting or downloading prefabri- 
cated opacity tables and using them in environments where grain condensation is significant. 
It is important to include the physical effects that are desired to have accurate opacities, 
most notably changes in the size distribution can greatly affect the Rosseland mean. 

Opacity tables for the F05 base set are available for download at our web site at 



http://webs.wichita.edu/physics/opacity, Custom opacity tables for non-standard compo- 



sitions or for different physical assumptions are available upon request, as are tables in 
gas density or gas pressure space in stead of log R. Note that a full set of opacity tables 



as outlined in iFerguson et all (120051 ) with (now) 155 values of hydrogen, X, and metal, Z, 
abundances take approximately a week of computer time, due to our full treatment of the 
molecular opacity sampling. 
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